A formula for charmonium suppression 
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In this work a formula for charmonium suppression obtained by Matsui in 1989 
is analytically generalized for the case of complex cc potential described by a 3- 
dimensional and isotropic time-dependent harmonic oscillator (THO). It is suggested 
that under certain conditions the formula can be applied to describe J ftp suppres- 
sion in heavy-ion collisions at CERN-SPS, RHIC, and LHC with the advantage of 
analytical tractability. 



I. INTRODUCTION 



The modification of the charmonium production cross section has been studied using a 
schematic 3-dimensional harmonic oscillator for the intermediate and final cc pair in 
In that reference the distorted wave Born approximation was used for the two-gluon fu- 
sion model and suppression ratios were calculated. In the present paper, we consider a 
3-dimensional THO with a complex and continuous time dependent frequency. For such a 
generalization, we derive the suppression ratio for charmonia states and present a formula 
for J/tp suppression including feed-down contributions. 



II. QUANTUM MECHANICAL EVOLUTION OF THE cc STATE 

The Charmonium suppression ratio was denned as a ratio of two cross sections by the 
expression S^(t) = ^^g^yj) anc ^ was calculated explicitly in Ref. jlj. From Eqs. (2.22) and 
(4.17) of that paper the survival probability for the s-wave can be written in the following 
form 
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SJt) 



J °° dr r 2 ip(r) U C c(r,t) 



lim Jq 00 dr r 2 ip(r) U C c(r, t) 



(1) 



III. TIME EVOLUTION OPERATOR FOR THE THO MODEL 



We make use of the standard path integral approach in order to calculate the time evolu- 
tion operator U c5 (r,t). We start by considering a 3-dimensional isotropic THO model with 

+ # w 2 (r) r 2 (r) , where r is the cc separation and the complex 



the Hamiltonian H 
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function of time uj(t) enters in the classical equation of motion for the heavy pair as 

r(r) + uj 2 (t) r(r) = . 



(2) 



The general solution of equation (J2J) is a linear combination given by r(r) = 

p(r)(^A cos7(r) + B sin 7(7-)), where 7(7-) = Jq dt' ^Ty- Replacing these definitions into 
(T5]), clearly leads to the following Ermakov equation |2] 

1 



p(r) + u (r) p(-, 



p\r) 



. 



If t G [0 , t] then A and B can be easily obtained from the initial conditions as 



.4 



r(0) 



B 



r(t) r(0) 

W)~W) 



cos7(i) 



(3) 



(4) 



p(0) ' sin 7 (t) 
Where we have used that 7(0) = 0. By replacing A and B in the general solution, we 
obtain r(r) and r(r). For a THO the classical action s c i and the fluctuation factor F(t) in 
the 3-dimensional isotropic space are defined in Ref. jsj. We calculate here their relationship 
with Ermakov function as 1 



d = £ (r(t) f(t) - r(0) r(0)) 



2 
/£ 

2 shi7(t) 



x 



r(t) 2 ( 7 (t)cos 7 (*) + ^| sin 7W 



+r(0) 2 (7(0) cos7(t) 



P(0) 
P(0) 



Pit) 

sin 7 (t))-r(t)r(0) 7(f) + 



pit) 



dr(t) 



JL( 



p pit) at) 



\ 2-ki p(0) sin 7(f) 



(5) 
(6) 



1 We use the notation r(t) — dr ^ \ T=t for all functions of time. 
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The time evolution operator for THO is given exactly by U(r,t) = F(t) exp(i s c /). In 
the present context, it will represent the quantum mechanical evolution of a cc state for a 
medium-modified (distorted) interaction up to the time t when it gets projected onto the 
asymptotic bound state spectrum. Thus we define U C c(r,t) = U(r,t). In fact formula (CQ) is 
independent of the initial condition which may be taken as r(0) = 0. 



IV. THE THO FORMULA FOR CHARMONIUM SUPPRESSION 



The ground state of charmonium J / ip can be identified with the Is- wave of the harmonic 
oscillator given by ip(r) = ijj{0) exp(-=^-j with 
gaussian shape over r appearing in 



4j. Thus we integrate the 



■ / V & U * 

which leads to the following suppression 



P(t) 



P(P) 



cos7(t) + 



/p(t)p(t) 1 



+ ^ 



(7) 



V j(t) j(t) 

The formula ((Tj) depends on j(t), the frequency and the Ermakov function 
p(t). For the case of the charmonium state ip' we take the 2s- wave given by <p(r) = 
| tp(0) — exp(;=^-J. Applying the formula jTJ we obtain 



Sf(t) — Sj/^(t) 



2 i sin 7(2) 



(8) 



For the Charmonium state Xc we take the 2p-wave given by x( r ) = x'(0) r expf^^-j. 
However, in this case there is a contribution of the angular momentum and it was shown in 
Ref. |l| that for such waves the formula ([TJ vanishes and the next-to-leading order term in 
momentum 0(p/m) must be considered leading to the expression 



s x (t) 



f °° drr* X (r) U' c5 (r,t) 



(9) 



with U> 5 



fi r 



2sm 7 (t)'(' ( t ) 7(0 P( Q ) X +P(°) 7(0) pit) x jUcc . The observable J/ip suppression 
ratio is influenced by feed-down from the higher charmonia states and we shall assume the 
following composition of the total contribution 



S{t) = 0.6 Sj/^t) + 0.3 S x {t) + 0.1 S#(t) . 



(10) 



The case of no feed-down is described by the expression S no (t) = Sj/^it) . Since we 
have already shown that S x (t) < Sj/^(t) and S^(t) < Sji^it) for Sj/^(t) < 1 it is clear that 
S(t) < SUt). 
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V. SUMMARY 

We have generalized Matsui's harmonic oscillator model for charmonium suppression to 
the case of time-dependent complex oscillator strengths and included the effects of feed-down 



on the J/ip suppression ratio. Preliminary results 
results from CERN SPS and RHIC can be found in 



:or the comparison with experimental 

3. 
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